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Abstract. We study filling phenomena in 3D wedge geometries paying 
particular attention to the role played by a line tension associated with the 
wedge bottom. Our study is based on transfer matrix analysis of an effective one 
dimensional model of 3D filling which accounts for the breather-mode excitations 
of the interfacial height. The transition may be first-order or continuous (critical) 
depending on the strength of the line tension associated with the wedge bottom. 
Exact results are reported for the interfacial properties near filling with both 
short-ranged (contact) forces and also van der Waals interactions. For sufficiently 
short-ranged forces we show the lines of critical and first-order filling meet at 
a tricritical point. This contrasts with the case of dispersion forces for which 
the lines meet at a critical end-point. Our transfer matrix analysis is compared 
with generalized random-walk arguments based on a necklace model and is shown 
to be a thermodynamically consistent description of fluctuation effects at filling. 
Connections with the predictions of conformal invariance for droplet shapes in 
wedges is also made. 
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1. Introduction 

Fluid adsorption on micropatterned and sculpted solid substrates exhibit novel phase 
transitions compared to wetting behaviour at planar, homogeneous walls [H El E]. 
The simple 3D wedge geometry has been extensively studied in the past decade 
theoretically [l[5lll[Zl[i[3[10l[lTl[ia[l3l[nj, experimentally [15l S [M E] and by 
computer simulation [18l [191 HH] • Thermodynamic arguments [211 [111 [23] show that 
the wedge is completely filled with liquid provided the contact angle is less than 
the tilt angle a. These studies show that the conditions for continuous wedge filling 
transition are less restrictive than for critical wetting at planar walls [SI [5] . Close to 
critical filling, the substrate geometry enhances interfacial fluctuations, which become 
highly anisotropic. We refer to these as breather modes excitations |5l[6]. However, 
most of these studies neglect the presence of a line tension associated with the wedge 
bottom. Previous studies by the authors [131 E] for short-ranged binding potentials 
show that the line tension may play an important role in filling phenomena and may 
drive the transition first-order if it exceeds a threshold value. We extend our analysis 
to arbitrary binding potentials, in particular to van der Waals dispersive interactions. 
Again, this shows that we can induce first-order filling by tailoring (micro-patterning) 
the substrate close to the wedge bottom. This may provide a practical means of 
reducing the fiuctuation effects which would otherwise dominate any continuous filling 
transition. This finding may have technological implications for microfluidic devices. 
However, the borderline between first-order and critical filling depends on the specific 
range of the interactions. If the binding potential between the interface and the 
flat wall decays faster than l/z"*, where z is the local interfacial height above the 
substrate, both regimes are separated by a tricritical point, as in the case of contact 
binding potentials [131ll4j. On the other hand, for longer-ranged binding potentials, a 
critical end point separates the flrst-order and critical flUing transitions. We note that 
these two situations correspond exactly to the fluctuation-dominated and mean-field 
regimes for critical filling [5] W ■ 

Our Paper is arranged as follows: In Section II we review bricfiy the 
phenomenology of wedge filling and introduce the breather mode interfacial model 
used in our study. The definition of the path integral used in our transfer matrix 
analysis is discussed in some detail. While other formalisms have been forwarded they 
suffer from a number of problems. As we shall show our definition is consistent with 
thermodynamic requirements (exact sum-rules), generalized random walk arguments 
and also the predictions of conformal invariance. Section HI is devoted to the analysis 
of wedge filling for contact binding potentials. Some of these results have been 
previously reported, without derivation, in a brief communication |13l I14j. Section 
IV extends the transfer matrix analysis to the important practical case of filling with 
long-ranged van der Waals forces. We conclude with a brief discussion and summary. 

2. The model 

Our starting point is the interfacial Hamiltonian pertinent to filling in shallow wedges 
(small tilt angle a) [4J: 

^{Vzf + Wiz-a\xlx) (1) 
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Figure 1. Schematic illustration of a typical interfacial configuration in the 
3D wedge geometry. The diverging lengthscales at the filling transition are 
highlighted. 



where z{x, y) is the local height of the liquid- vapour interface relative to the horizontal, 
E is the liquid-vapour surface tension and W{z, x) is the binding potential between 
the liquid- vapour interface and the substrate (see figure [1]). Here after we assume the 
temperature defines the energy scale and set fc^T = 1. To first approximation we 
may suppose that W{z, x) is independent of the position across the wedge section, i.e. 
W{z, x) K, Wt^{z), where Wt^{z) is the binding potential due to a single flat substrate. 
Corrections may arise when the liquid adsorption is small enough however - a point 
we shall return to later. 

A mean-field analysis shows that locally the interface across the wedge is fiat and 
that fluctuation effects are dominated by pseudo-one-dimensional local translations in 
the height of the flUed region along the wedge (the breather modes) [SJIS]. Fluctuation 
effects at filling can be studied by using an effective pseudo-one-dimensional wedge 
Hamiltonian which accounts only for the breather-mode excitations [Sj [6] : 

nw[l]^ I dy[^{^\\vw{l)\ (2) 



dy 

where l{y) = z{0,y) is the local height of the interface above the wedge bottom. The 
effective bending term A{1) resisting fluctuations along the wedge can be expressed as 

urn-. 

2^1 2t 2U 

A(/) « + — « (3) 

a a a 

where r is the line tension associated with the contact lines between the filled region 
and the substrate far from the wedge bottom. Note that for large I, we may neglect 
the r contribution since A(^) is proportional to the local interfacial height. Similarly, 
for large I the effective binding potential Vw{l) takes the form |5ii6j: 

a a 
+2t + t'+ dxWnil-a\x\) (4) 

The first two terms corresponds to the bulk and surface thermodynamic contributions 
required to form the filled liquid region [7]. Here h denotes the bulk ordering field 
measuring deviations from bulk two-phase coexistence, 9 is the contact angle of the 
liquid drop at the planar wall-vapour interface and Zjr is the equilibrium liquid layer 
thickness for a single planar wall. The line tension r is defined as above, and t' is the 
line tension associated with the wedge bottom. Note that the line tension contributions 



Controlling the order of wedge filling transitions: the role of line tension 



4 



are essentially independent of I for ^ ^ Z^, so they become irrelevant in that limit. 
Finally, the last term corresponds to the binding potential contribution to Vw- Upon 
minimisation of Vw{l), we recover the mean- field expression for the mid-point height 
(at bulk coexistence) [SI : 



V^2 

W^il) + — = AW^il) 



(5) 



However, we stress that the form (j4]) for Vw{l) is only valid for I l^. For I ^ Itt, both 
A(/) and Vw{l) will behave in a different manner. Furthermore, we may control the 
adsorption properties for small / by micropatterning a stripe along the wedge bottom, 
so as to weaken the local wall-fluid intcrmolccular potential. The interfacial binding 
potential is consequently strengthened, and under some conditions it may bind the 
liquid-vapour interface to the wedge bottom even at the filling transition boundary 
9 = a. Thus by introducing a line tension associated with the wedge bottom one may 
induce first-order filling in the modified wedge provided the modification is strong 
enough. As we shall see the lines of first-order and continuous filling transitions are 
separated by either a tricritical point or a critical end point depending on the range 
of the intermolecular forces. 

The quasi-one-dimensional character of the Hamiltonian means it is amenable to 
a transfer-matrix analysis. The partition function corresponding to this Hamiltonian 
can be expressed as the following path integral [24] : 



Z{hJa,Y) = JvicM-^w[l]) 



(6) 



However, the presence of a position-dependent stiffness coefficient makes the definition 
of the partition function ambiguous. This problem was pointed out, but not 
satisfactorily resolved, in \B[ and is intimately related to issues associated with the 
canonical quantization of classical systems with a position-dependent mass [251 126j . 
In this paper we use the following definition of the partition function: 

N 

ZihJa,Y)= lim dh...dlN-iT\Kil,J,^uY/N) (7) 
where Iq = 1^ and In = li,, and K{1, 1', y) is defined as: 

The partition function la, Y) satisfies the differential equation 

dz{hx,Y) 



HwZ{lb, la, Y) 



dY 



(8) 



(9) 



with initial condition Z{li,, la, Y) — > d{li, — la) as 1" — > 0. The operator Hyy is defined 
as (Ml: 



where Vw{l) is given by 



ld_ 



_2 d_ 

A(4)% 



Vw{l) = - 



2A(0 



3 fA'{l) 



Vw{lb) + Vw{lb) 



A"(0 



4 VA(0 



2A(0 



(10) 



(11) 
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and prime denotes differentiation with respect to argument. Tlie solution of ^ can 
be expressed via the spectral expansion |25 ■ 

Z{h,la,Y) = ^Va(W:Ga)e-^°'' (12) 
a 

where {ipail)} is a complete orthonormal set of eigenfunctions of the Hamiltonian 
operator Hw, with associated eigenvalues Ea- 

Analogous to discussion of 2D wetting [21] we can now obtain the interfacial 
properties from the knowledge of the propagator Z(lb,la,Y). In particular, the 
probability distribution function (PDF), Pw{l), can be obtained as: 

while the joint probability Ppj/ (/i, Z2, ^l, ^^2) of finding the interface at midpoint 
heights li and I2 at positions Yi and Y2{> Yi), respectively, is: 

P'^\ll,l2,Yi,Y2) ^ Z{l2,luY2-Y^) 

ZihJ-L/2,Y,+L/2)Z{lL/2,l2,L/2~Y2) ^^^^ 

X lim ; — (14) 

1-^°° Z{Il/2J^L/2, L) 

Further simplifications arises if we assume the existence of a bounded ground 
eigenstate ipo. Then, substitution of the spectral expansion (I12p into (fT3|) reads (for 
an infinitely long wedge) 

Pw{l) = \Ml)\^ (15) 

Similarly the excess wedge free-energy per unit length of an infinitely long wedge 
is identified with the ground eigenvalue Eq. The two-point correlation function 
h{li, I2, Y) can be obtained in a similar manner: 

h{h,l2, Y) = P^'> {h,l2, 0, Y) - Pw {h)Pw ih) 

= ^a{h)Mh)Ml2Wo{l2) exp [^{E^ - E„)Y] (16) 

For large separations this vanishes exponentially allowing us to determine the 
correlation length = {Ei — Eo)~^, where Ei is eigenvalue corresponding to the 
first excited eigenstate. This result still holds even if Ei corresponds to the lower 
limit of the continuous part of the spectrum of Hw, although now the leading order 
of h{li, I2, Y) is not purely exponential but is modulated by a power of Y. 

An interesting connection with 2D wetting problems can now be seen. Introducing 
the change of variables [57] 

V{l)^jdl^) ; ^(77) = A(0-^/V (^(^)) (17) 

the eigenvalue problem Hw'ipai.l) — Eaipail) transforms to the following Schrodinger- 
like equation: 



2 dif 



+ {Vw[im + V^iv)) Mv) = EMv) (18) 



where is defined as: 



^wKV) 32A2(/(77)) \ dri ) 8A(?(r/)) d??2 ^ 
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(|18p shows that the fiUing problem can be mapped onto a 2D wetting problem in the rj 
variable, under an effective binding potential. In our case A{1) = 21^1 /a, so we obtain 
Vwil) = -3a/16E/3. The change of variables HID leads to ry = ^SS/Oa^^/^ and 
thus ~ — 5/72?7^. Note that the use of the variable 77 cx l^^^ as the appropriate 

collective coordinate has been previously recognized in the literature [51 [7]. However, 
our definition of the path integral ([7]) and ([H) leads to a novel term in the wedge 
binding potential, which will be essential in our study and ensures thermodynamic 
consistency. 

The filling potential Vw{l) must fulfill some requirements. In order that the 
interface does not penetrate the substrate, Vwil) has a hard- wall repulsion for 
^ < 0. Consequently we need to impose an appropriate boundary condition on the 
eigenfunctions of Hw at Z = 0. The analytical expression of the boundary condition 
is obtained by a regularization procedure: we assume that the filling potential Vw 
and the position-dependent stiffness A are constant for / < ^Oj where ^0 is some 
microscopic scale. Furthermore, we impose that A(Z) must be continuous at Z = ^0, so 
^{l < (,q) — 2S^o/"- On the other hand, the filling potential can take an arbitrary 
value —U. The latter square- well potential models the modification of the filling 
potential Vw{l) for small I due to the line tension associated with the wedge bottom. 
The eigenstates must fulfill the usual matching conditions that (j)a and d<j)a/dl are 
continuous at Z = ^o- Finally, we consider the appropriate scaling limit as ^0 ^ 0. 

The qualitative form of the filling potential depends on the order of the mean-field 
phase transition [5, 6 . For critical filling with long-ranged intermolecular forces, the 
binding potential at bulk coexistence behaves as 

vv(o^S?!^,^,-,,„^ ,20, 

a (p — l)a 

where ^ is a Hamaker constant while the exponent p depends on the range of the 
forces. Specifically, for non-retarded van der Waals forces p = 2. For systems with 
short-ranged forces this is replaced by an exponential decay ~ exp(— kZ) where k is 
an inverse bulk correlation length. The presence of the fluctuation-induced filling 
potential Vw{l) oc l~^ in ([9]) gives rise to two distinct scenarios. For p > A and 
large I, the direct contribution in ([5]) cx /^^^ is negligible compared to the fluctuation- 
induced potential Vw arising from the position dependent stiffness. In this case we 
anticipate universal, fluctuation dominated behaviour. On the other hand, for p < A 
and large I we can neglect Vw ■ Since the Z^"^ contribution to the binding potential 
is now repulsive, we expect a qualitatively different phase diagram. It is worthwhile 
noting that both situations exactly correspond to the existence of a mean-field and a 
filling fluctuation-dominated regime for the filling transition predicted from heuristic 
arguments [3 [6]. 

Fig. [2] shows the schematic filling phase diagrams we expect for short-ranged 
(Fig. 2a) and long-ranged forces (Fig. 2b). Previous work for the contact potential 
case [m [14] showed a similarity between 3D wedge filling (for h — Q) and 2D wetting, 
with 9 — a playing the role of the ordering field and the effective line tension the role of 
the wetting binding potential. The borderline between the first-order and the second- 
order transition lines corresponds to a tricritical point (analogous to the 2D critical 
wetting case). For long-ranged forces the analogy between 3D filling and 2D wetting 
still holds (see above). However, as the effective binding potential is repulsive for 
large I aX 6 = a, the similarity must be established with 2D first-order wetting. As for 
contact binding potentials, we expect the interface to be bound to the wedge bottom 
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Figure 2. Schematic phase diagrams for 3D wedge fiUing pertinent to: (a) short- 
range forces, (b) long-range forces. U is the effective hne tension strength, 9 is the 
contact angle for the liquid on a planar substrate, and a the wedge angle. The 
filled circle locates the borderline between the first-order transition line (dashed 
line) and the second-order transition line (thick continuous line) between the 
bound and unbound states. These correspond to a tricritical point and a critical 
end-point for short-ranged and long-ranged forces, respectively. In the latter case, 
the first-order transition line continues to a first-order pseudo-transition line in 
the partial filling region (dotted line), which terminates at a pseudo-critical point 
(cross). See text for explanation. 



at 6 ~ a for a large well depth U (corresponding to the line-tension contribution) . As 
U decreases, the interface will unbind along the 9 = a path. However, the interface 
must tunnel through a free-energy barrier to become unbound. Consequently, the 
borderline between the first-order and second-order wedge filling is a critical end-point, 
where the spectator phase is the bound state at 9 — a. Actually, the connection with 
wetting phenomena also leads naturally to this picture, since first-order wetting was 
previously recognized as an interfacial critical-end point scenario |28j . In principle, we 
may expect the first-order transition line to continue in the partially filled region as 
the coexistence between two bound states with different adsorptions. This thin-thick 
transition is analogue to the prewetting line, and it should terminate at a critical 
point. However, the quasi-one-dimensional character of the wedge geometry rules out 
this transition, as it is destroyed by breather-mode fluctuations. Nevertheless traces 
of this smeared transition may be found in the bimodal form of the interfacial height 
PDF (see later). 

We will consider two cases which correspond to different regimes of the filling 
transition. The contact interaction will be studied as the paradigm of the filling 
fluctuation-dominated regime. On the other hand, the van der Waals (p — 2) case is 
analysed as a prototypical case of the mean-field regime. 

3. Results for contact interactions 

We consider first the case of short-ranged potentials. Some of our results have been 
reported in a brief communication but without detailed explanation [131 114| . Here we 
provide full details of our transfer matrix solution and present new results for the form 
of the propagator. At lengthscales much larger than that of the bulk correlation length 
we may write Vw{l) = S(6'^ — a'^)l/a for Z > and allow for line tension arising 
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Figure 3. Scaled probability distribution function Pwi^/iu) as a function of 
the scaled wedge midpoint interfacial height l/^u for 6 = a. Inset: Plot of 
gi/3Ai'(ei/3)/Ai(e:-'-/'^) (continuous line) and \/u — ecot — e — 1/2 for u = 
0.5 < Uc (dotted line), u = Uc = 1.358 (dashed line) and u = 2.0 > «c (dot- 
dashed line). For the latter case, the value of the reduced ground eigenvalue 
e = e* is highlighted. 



from the wedge bottom via a suitable boundary condition at the origin. Analysis of 
^ shows that the short-distance behaviour of the eigenfunctions will be dominated 
by the l^^ contribution to the effective filling potential. In fact, ipail) ^ l^^^ or l^^^ 
as Z — > 0, so the PDF function Piy{1,9) ~ Z or as Z — > 0. We anticipate that the 
latter behaviour corresponds to the critical filling transition as predicted by scaling 
arguments [7], but the former will be a completely new situation which, as we shall see, 
is related to the possibility of tricriticality. Turning to the critical filling transition 
first we note that the short-distance behaviour of the PDF, Pw{l,0) ~ P, which 
emerges from our analysis, is indeed the required, thermodynamically consistent, 
result. This short-distance expansion ensures that the local density of matter near 
the wedge bottom contains a scaling contribution that vanishes cx T — Tj, where Tf is 
the filling temperature. This is the required singularity which emerges from an analysis 
of sum-rules connecting the local density near the wedge bottom to (derivatives of) 
the excess wedge free-energy [29 . This leads us to conclude that our definition of the 
path integral in ([5]) is the correct one for the 3D wedge filling problem. As we shall 
see it also ensures that our model is conformally invariant. 

3.1. Wedge filling along the 9 = a path 

At the filling phase boundary {6 — a), the problem can be mapped onto the 
intermediate fluctuation regime of 2D wetting [30| via (fT8|) . We apply the 
regularization method described above, and consider the scaling limit — > 0, L/ — > oo 
and 4EC/^Q/a u. Under these conditions, the ground eigenvalue Eq satisfies 

, , f w,Ai'(e~V3) 

V^cot - - = '^"'^4^ (21) 
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where e = —AYiEq^^q /a and Ai(a;) is the Airy function. Graphical solution of (PTjl 
(see inset of figure ^ shows that there is a bound state with i?o < for u > Uc, 
where Uc ~ 1.358. Otherwise Eq = and there is no bound state. The existence 
of a bounded ground state a,t — a implies that the filling transition is first-order 
for u > Uc, and critical for u < Uc- The explicit form of the PDF for u > Uc in the 
thermodynamic limit is 

6%/37r I 



Pw{l,9^ a) 



(22) 



Note that the lengthscale ^„ can be 



where 

arbitrary large as u ^ Uc- Figure [3] plots the PDF in terms of the scaling variable 
l/^u- For small I, Pw{l,ct) ~ P", with a short-distance exponent (SDE) ju = 1. 
Asymptotically Pw{l,a) ~ i/Iexp [-4(Z/^„)^/^/3] as Z ^ oo. 

The mean interfacial mid-point height Iw = {1} and roughness = ^ (P) — (l)"^ 
satisfy Iw ^ £.± ^ showing that, in the scaling limit, there is only one lengthscale 
controlling the fluctuations of the interfacial height. Finally, the correlation length 
along the wedge axis close to the filling transition can be obtained as = — l/£^o = 
4E^^/ck cx (m — Uc)~^ since we can identify Ei = 0. 

These observations indicate the emergence of a new relevant field (in the 
renormalization group sense) tu oa {u — Uc), in addition to tg oc 6 — a and the bulk 
ordering field h. Thus the conditions 9 = a, u = Uc and h = correspond to a 
tricritical point which separates the lines of first-order and critical filling transitions. 
The excess wedge free energy density Eq vanishes as Eq ^ tu ~ t„ as u tends to 
Uc from above. Critical exponents for the divergence of the characteristic lengthscales 
can be defined as the tricritical point is approached along the 9 = a path: 



I'll 



tn 



6 



tu 



til 



(23) 



Our results show that = i^" = 1 and i/^ = 3. Finally, the effective wedge wandering 
exponent ^vv — i'±/'^y = V^: which coincides with its value for the critical filling 
transition [S1[S]. 

For u < Uc and 9 — a, the interface is unbounded in the thermodynamic limit. 
However, we can study the finite-size behaviour of the droplet shape when the interface 
is pinned very close to the wedge bottom at positions y — -izL/i. Making use of results 
presented in the Appendix, in particular (IA.11|) and (|A.12p . we find that the PDF at 
the tricritical (u = Un 



and critical (u ^ Uc) wedge filling transition are given by 



3i/3r(2/3) 



exp 



pW^Uc 



\l,9^a) 



Xll^ 



35/3r(4/3) 



exp 



9 

(Al/)^ 



(24) 



where 



Xr 



( aL ) 



1/3 



1 



2Y 
T 



-1/3 



(25) 



(26) 



The typical droplet shape may be characterized by the most-probable position lmp{y) 
which follows from the relation dP^{lmp)/dl = 0, or by the average shape lav{y) via 
the definition = dllP^{l) [24]. In all cases we find that the typical droplet 
shape follows XlI = c, where c is a number which depends on the definition of the 
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typical shape and whether the pinning is at critical or tricritical filling. Consequently, 
the droplet shape must obey: 

Y^ + ^l^^^ (27) 

Satisfyingly this shape is precisely that predicted by the requirement of conformal 
invariance which can be used to map a droplet pinned at just one end to one pinned 
at two points {y = ±L/2) [21]. This gives further indication that the definition of the 
measure used in our transfer matrix formulation is appropriate for the wedge geometry. 

3.2. The necklace model 

A simple model can be introduced to understand the filling properties of the wedge at 
9 = a. This model is a generalization of the necklace model introduced for 2D wetting 
[SH [33] . The interface is pinned to the wedge bottom along segments of varying length 
but unbinds between them, forming liquid droplets. We associate with each vertex a 
weight V which is related to the point tension between a bound and unbound state. 
Following the analysis described in [32] , we introduce the generating function: 

oc 

G(z) - -^^^i (28) 

where Zj^ is the interfacial canonical partition function of a segment of length La 
(which we assume to be discretized in intervals of length a) and z is an activity-like 
variable. The pure pinned (A) and unbound {B) states have generating functions: 

G^ = ^Zfz^~^(.z)^ (29) 

L L 

Gb=E^^"-'?"E^ (30) 

L L 

where s — exp(— u), with u as an effective line tension in units of ksT/a. On the 
other hand, w = exp(— tq), where tq is the reduced excess free energy per unit length 
of a long liquid droplet on the wedge, which we can assume to be zero by shifting the 
energy origin. Finally, "0 is the exponent characterizing the first return of the interface 
to the wedge bottom. This can be calculated from the u Uc limit corresponding 
to the completely filled regime. As shown in (jA.13|) and (|A.15[) in the Appendix, 
■0 = 4/3. Now the complete generating function is: 

G{z) = Ga + GavGbvGa + ... = , , (31) 

1 - v^Ga[z)Gb{z) 

Thus there will be a continuous phase transition at Uc defined as [32] : 

uc = -lii[w{l-v'^Gc)] (32) 

where Gc — /^^ ■ Below Uc, the interface is completely unbound while it 

remains pinned for u > Uc- The singularities of the various interfacial properties 
close to Uc are characterized by critical exponents. In particular, the specific heat 
critical exponent and the longitudinal correlation length critical exponent Vy can 
be expressed in terms of the exponent -0 as [311 133] 

2-alV=^" = ^ (33) 
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Figure 4. Plot of the scaled PDF for to = —1.5 (thick dashed line), eo ~ 1-086, 
which corresponds to u = Uc (thick continuous line) and eo ~ 1.639, corresponding 
to M <^ Mc (thick dot-dashed line). For comparison, the PDF from (1221 with 
~ 1. 968^6 (which corresponds to eg = —1.5, see inset) is also plotted (thin 
dashed line). Finally, the scaled PDF obtained in [5] is also shown (thin dot- 
dashed line). Inset: Plot of the ground (black lines) and first-excited (red lines) 
reduced eigenvalues e = TiE^g/a of Hur as a function of for u < Uc 

(continuous lines) and u > (dashed lines). 



Substituting i/; = 4/3 we find the correct values 2 — = — 3 derived earlier. In 
order to estimate the interfacial average height and roughness, we need the interfacial 
height PDF in a hquid bubble, which is given by pS]). A similar argument to the one 
presented in [32] leads to: 



^4/3 ^7/3 

, (34) 

where Yb is the length of a liquid bubble, and a is the average of a in the pure liquid (B) 
phase. It is straightforward to see that Iw ^ ^ fy^*^, with (w — 1/3. Consequently, 
^± — l^w = I7 in agreement with our exact results. Note that in contrast to the 
2D wetting case [32l [33] the exponent associated to the probability of first return 
tfj 2 — (w- This can be traced to the role played by the position dependent stiffness 
in the filling model ^ which biases the random-walk-like motion of the interface. As 
pointed out previously [14] it is remarkable that the critical exponents for 3D critical 
and tricritical wedge filling are identical to those anticipated for 2D complete and 
critical wetting with random bond disorder. The necklace model provides an elegant 
means of understanding this unusual dimensional reduction. 

3. 3. Wedge filling for 9 > a 

Now we extend some of our previous results to partial filling conditions, i.e. 6 > a. The 
rescaling = (E/a)i/4(6»2 - a^)^/^ri/V2, = {Y./a)-^/^{9^ - a^)-^/*Ea transforms 
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the Schrodinger equation (|18p into a parameter free form 



4 dry2 



3~x2/3 



144^2 



(35) 



This leads directly to the critical behaviour of the mean mid-point height Iw = [I) ^ 
{6 — a)''^/^ and the correlation length along the wedge iy ^ {d — a)"'^/'*, in agreement 
with scaling arguments [SI [5], provided that the PDF is not singular, i.e. takes non- 
negligible values at finite values of r). The PDF is obtained analogous to the case 
= a case described above, from determination of the ground eigenstate We find 
for the PDF 

/ 



Pw{l^ > a) — C-^ exp 



4 



where 



H,2 

4 2 

S-i/2[(0/a)2 



I 



11-1/4 , — 

i , en = 



(36) 



S-Eo'^l/a is the reduced ground eigenvalue, 
C is a normalization factor and Hs{z) is the s-order Hermite function [34]. As 
I — > 0, the PDF vanishes (in general) like Pw{l,0) ^ I while at large distances 
Pwil,0) - ?'o/2-iexp[2?(eo - l/^9)/^e] as / ^ oo. 

The value of the reduced ground eigenvalue eg depends on the boundary condition 
at I = 0. In order to do this consistently we once again turn to a regularization 
procedure: in the scaling hmit (fo and 4I][/^Q/a —> u), the reduced eigenvalues 
are the solutions of the equation: 



/u cotvu - 



V2 



V2 



V2 



(37) 



where eg is the minimum value of them. For u close to Uc, we can expand the left-hand 
side of ([57)1 around Uc as: 



?0 



/u COtVM — 



(38) 



where the positive (negative) sign corresponds to u > Uc (u < Uc), respectively. This 
expression allows us to identify ^„ oc S,o/\u — Uc\ in a manner completely consistent 
with our expression obtained for u > Uc &t = a. 

We obtain scaling behaviour for the wedge excess free-energy per unit length i?o 
(see also inset of figure H]): 



Eo 



{9 -a) 



3/4 ± 
^0 



■(6'-a)i/4 



(39) 



where c is an unimportant metric factor, and the sign corresponds to the situations 
u > Uc and u < Uc as above. The scaling two functions and €q have the following 
properties: 



(0) 



e+(x- 



6(7(0) = 1.086 (40) 
+oo) - -x^ (41) 
(x +c») = 1.639 (42) 

The asymptotic behaviour of the PDF as the filling transition is approached, i.e. 
> a, is different for three situations: (i) u > Uc, (ii) u < Uc and (iii) u = Uc 
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(see figure S]). For the case (i), saddle-point asymptotic techniques [3S] applied to 
the PDF recover the expression for the PDF (|^^ . Consequently, the lengthscale 
governing the interfacial height and range of the breather-mode fluctuations is 
which remains finite as S,d diverges. Thus on lengthscales compared to the PDF 
becomes a highly localized delta function located at / = 0. On the other hand, the first- 
excited eigenvalue scales as Ei cx (see inset of figure HJ, so the lateral correlation 
length also remains finite, oc 

For case (ii) we must take the limit (j42|) of the scaling function Cq . Substitution 
into (|36p leads to the asymptotic behaviour of the PDF as 9 ^ a. Now, the 
lengthscale which controls both the mean interfacial height and roughness is ^e- It 
is also interesting to note that Pw{l) ~ for / 0, so thermodynamic consistency 
is assured. Our solution is different from the PDF reported in [H 1^ (see figure H]) 
although the global behaviour is qualitatively similar. On the other hand, the lateral 
correlation length has the asymptotic behaviour £_y = {Ei — Eo)~^ ^ £,g. 

Finally, the PDF for case (iii) is obtained by substitution of the limiting value 
(|40p for Co into (|36p. Although the relevant lengthscales behave asymptotically as in 
the case (ii), the scaled PDF is different, as shown in figure[4l In particular, Pw{l) ~ I 
as Z ^ 0. 

We can define the critical exponents for fixed u as: 

Eo ^ t]-""- , l^ ^ i/"- , e± ^ t-/^ , - tg"' (43) 

where we define tg cx 9 — a. For both critical and tricritical filling, our analysis 
shows that the critical exponents take the values 2 — a^, =3/4, ~ = 1/4 
and Vy = 3/4, so the wandering exponent is Cw = '^±/'^y = V^- These findings are 
in agreement with scaling predictions 5,f6]. Finally, we note that the tricritical gap 
exponent A*, which relates the critical exponents along the 9 = a and fixed-u paths, 
can be obtained from the scaling form of Eq (|39|1 as A* = 4 [14] . 



4. Results for dispersion forces 

Wedge filling in systems with short-ranged forces is representative of the universality 
class of fluctuation dominated behaviour occurring if the exponent in the binding 
potential p > A. In almost all practical realizations of wedge filling however, 
long-ranged, van der Waals forces from the fluid-fluid and solid-fluid intermolecular 
potentials will be present. While exceptions to this may be found, for example in 
polymer systems, the case of long-ranged forces is much more the rule than the 
exception. We wish to understand how such long-ranged forces alter the nature 
of the phase diagram and in particular the change from first-order to continuous 
filling behaviour. We can allow for the presence of long-ranged forces through the 
binding potential (|20p . For the three-dimensional case and non-retarded van der 
Waals interactions the value of p = 2. While higher order-terms are present these 
will not play a significant role in determining the physics and can be safely ignored. 
Here we show that the potential 

.,, ^{9^~a^)l 2A 

Vwil) = — + -J (44) 

a at 

is amenable to exact analysis. The form of the total effective binding potential 
Vw{l) + Vwil) is shown in figure El As one can see there is a local maximum for 
small I which arises from the competition between the fiuctuation-induccd component 
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» 




l*=(E/a) 



Figure 5. The effective wedge binding potential Vw{l) + Vw{l) for an effective 



Hamaker constant 2 A/ a = 1 and: 
(dashed line) and (iii) 6^ — = - 



(i) 9^ — = 0.1 (continuous line), (ii) 
0.1 (dot-dashed line). 



and the van der Waals contribution. A rough estimate for its location can be obtained 
by setting 9 = a, and a simple calculation leads to 

r = (45) 

For 9 > a, a local minimum is obtained for larger I. This minimum arises from 
the balance between the thermodynamic contribution — 0^)1 /a and the van der 
Waals component of the effective binding potential. The position of the minimum is 
governed for — a <C a by the mean- field interfacial height: 



iMF _ 



2A 



Since for long-ranged forces we anticipate that mean-field theory describes correctly 
the critical filling transition, a third Icngthscale is the mean-field roughness, identified 
as [5l[6]: 



^1/4 



(47) 



(48) 



Since in our case Wj^ll) = —A/P, we find, from (l46 

where our definition of is unchanged from the previous Section. It is remarkable 
that the roughness is independent of the strength of van der Waals interactions. 

The fluctuation-induced contribution to the wedge binding potential suggests it 
may be possible to find a bound state at a lengthscale determined by I*, for 9 > a. 
However no bound state is possible for 9 < a: any interface would tunnel through the 
barrier and unbind completely from the wedge bottom. 

Our analysis proceeds along the same lines as the previous Section. First we 
investigate allowed states that exist at the filling phase boundary 9 = a, and then 
extend our analysis to the partial filling regime 9 > a. 
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Figure 6. Scaled probability distribution function Pw(l/l*) for 
u Uc in presence of dispersion forces. 



a and 



4-1- Wedge filling along the 9 — a path 

First we suppose 9 = a and search for a bound state V'olO of Hw- As anticipated this 
is only possible if E'o < 0, in which case the eigenfunction is given by 



(49) 



where we define as previously £,u = (— 4i;_Eo/ck)~^/'^. 

In order to obtain Eq, we make use of the regularization procedure. In the scaling 
limit this satisfies the equation: 



\/u — e cot \/u — e — - = - 



Ai' (ei/3 



9;; 

4(Z*) 



A 



/3 



9il 
4(FF 



(50) 



where e = {^o/^u)^- For values of e not too close to zero, ~ Co ^ I* in the scaling 
limit, so the second term in the Airy function argument can be neglected. Thus 
the solution for e coincides with that corresponding to a contact binding potential. 
However, as e — > 0, a crossover to a different situation is observed. In particular, for 
Cu ^ I* we can make use of the asymptotic expansion of the Airy function for large 
arguments: 



Ai(2;) 



exp(-2a;3/2/3) 



2^^x1/4 

Substituting into (|50|) . we find that 



+00 



(51) 



\/it — e cot Vi 



(52) 



3^0 II* f l* + 3Co 
21* 3^0 V 3Co 

Next we define Uc as the value of u at which e = 0. It is straightforward to see that 
Uc = 1.358 + 0(^0/^*), so in the scaling limit the threshold for the existence of a 
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bound state is the same as for the contact binding potential. Now we expand for 
u around Uf, and e around zero. In the scahng hmit, we obtain that: 

^~~(^("-"^) (53) 

Finally, the PDF in this regime can obtained from substitution of the asymptotic 
relationship ([5T|) into After some algebra, the PDF reads: 

^H'(0=(^cxp(-^) (54) 

(see also figure |6l). Since /* remains finite at Uc, the interface remains bound as u 
approaches Uc from above, and suddenly unbinds for u < Uc- This identifies the 
point h = 0, 6 = a and u = Uc as a critical end point in the surface phase diagram. 
However, although u is not a relevant field in the renormalization-group sense, the 
longitudinal correlation length = jii'ol"^ diverges as {u — Uc)~^. A similar behaviour 
occurs within the subregime C of the intermediate fluctuation regime for 2D wetting 
transitions [3D]. 



4-2. Wedge filling for > a 

We now extend our results to 9 > a proceeding in the same way we did for contact 
interactions. Again we can analytically obtain the ground state "00(0 as: 



ipoil) oc Vl 



exp 



I f I 



" 16(i*)^ 

where eo, and Hs{x) are defined as in p6)) for the contact binding potential. Note 
that the dependence on the dispersion forces only appears at the order s of the Hermite 
function Hs- 

The regularization procedure leads to the following equation for the reduced 
eigenvalues e (for u close to Uc)- 

' 'u cot V It w ±- 



?o V 2; r[i] ^„ 21* 

(-75) 



V2 



V2~ ) ^'^'^ ; _ (56) 

4 2 16(1. )2 ^ ' 

and eo is the minimum of the solutions. We can see that ep will depend now not only 
on the ratio Ce/^u and the sign of u — Uc, but also on £,0/1*. Nevertheless as for the 
case with a contact binding potential, we are mainly interested in the limit — > a, so 
that the lengthscale is very large. 

For u > Uc, saddle-point asymptotic techniques analogous to those applied in the 
previous Section [35] show that the ground state of Hw converges to the expression 
(|49p as — *■ a. Consequently we recover the results obtained earlier for the special 
case = a. This indicates that, for these values of m, the wedge filling transition must 
be first-order. 

For u < Uc the non-existence of ground state for 9 = a indicates that the filling 
transition is critical. We anticipate that mean-field theory will describe faithfully 
singularities at the critical wedge filling [S] [5] . Thus we expect that the interfacial 
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height PDF is centered around l^^ with Gaussian fluctuations on the scale of 
representing the breather mode excitations. The mean-field value of the excess free 
energy per unit length is given by Vwil^Y^)- Consequently, the mean-field value of 
the reduced ground eigenvalue Eq^^ is: 

e'o'' = § (57) 

The shift of eq with respect to eg^^ due to breather-mode fiuctuations Aeo can be 
estimated in the following way. We expand Vw around its minimum up to quadratic 
order, so the shift may be estimated via: 

Ai?o - iw^(^r)(er (58) 



implying that 



2 

~ 4f ^ ^ (59) 

which vanishes as 9 ^ a. 

We can now proceed with a more formal derivation. Equation (j55p can be written 

as: 

^^0(0 oc \/^exp(-a;V2)fl"s(x) (60) 
where we have defined: 

V2AI Aeo 



(61) 
(62) 

with Al = I - l^^ and Aeo = eo - e^^^. Note that is the wavefunction of the 
harmonic oscillator in the x coordinate (in units of yjh/muj), modulated by the factor 
^/l. As in the latter case, ■00 will increase exponentially asx— >— oo,i.e. Z^O and 
large l^^ ^ unless s is a non- negative integer. This result is independent of the explicit 
value of £,e/£,u- Consequently, the shift Ae for the lowest eigenvalues is given by: 



Ae = ^lie^r + 4 ( n + i) - - ^ (63) 



with n a non-negative integer. The ground eigenstate will correspond to the case 
n = 0. The corresponding PDF becomes a Gaussian: 



Pwil)^J^e^j>[-^^L^] (64) 



with (/) = /^^ + 2l*/3 and roughness = Finally, the longitudinal 

correlation length = {Ei - Eq)'^ ~ Ceeo^^/2 (6* - a)-^. 

Consequently, the explicit transfer matrix solution is in complete agreement 
with the predicted mean-field values for the critical exponents when p = 2 [Sj [6]: 
2-al = 1/2, = 1/2, = 1 and i^l = 1/4. 

Finally, we search for the existence of the thin-thick transition line. Our 
calculations show that there is no sharp phase transition for 6 > a. However, we 
observe that the interfacial PDF becomes bimodal for u < Uc and some range of 
values of > a. We can identify a first-order pseudo-transition line when the areas 
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l*/^,0.4 



0.2 - 



0.03 



0, 








5 



Figure 7. Interfacial height PDF at the first-order pseudo-transition for ^g/l* = 
4.216 (continuous line), 2.981 (dot-dashed line), 2.749 (dashed line) and 2.494 
(dotted line). Inset; location of the first-order pseudo-transition line. Note that 
ocUc — u (here u < u^) and (i*/^e)* oc 9 — a. The diamonds represent our 
calculated pseudo-coexistence values. The continuous line serves only as a guide 
for the eye. 



under each maximum of the PDF are equal (see Fig. [7|) . This hnc is the continuation 
to the partial filling region of the first-order filling transition line for u < Uc, and 
touches tangentially the filling transition borderline at u = Uc- As u decreases, the 
two maxima become closer, and eventually merge (the pseudo-critical point). 

5. Discussion and conclusions 

In this paper we have reported analytical results for 3D wedge filling transitions in the 
presence of short-ranged and long-ranged (van der Waals) interactions based on exact 
solution of the continuum transfer equations for a pseudo one-dimensional interfacial 
Hamiltonian. First-order and continuous (critical) filling are possible for both types of 
force depending on the strength of the line tension associated with the decoration of 
the wedge bottom. Our analytical solution for the interfacial height PDFs at critical 
filling recovers the known values of the critical exponents for critical filling for short 
and long-ranged forces. In addition we have elucidated the nature of the cross-over 
from first to critical filling which occurs at a specific value of the wedge bottom line 
tension. This is qualitatively different for systems with short and long ranged forces 
whose surface phase diagrams are shown to have tricritical and critical end-points 
respectively. To finish we argue that these results, which are proto-typical of the 
fluctuation dominated and mean-field regimes respectively, are qualitatively valid for 
any type of binding potential. As already mentioned, the 3D wedge filling phenomenon 
can be mapped onto a 2D wetting problem with a new collective coordinate rj cx Z^^^. 
If we suppose that, at the filling phase boundary, l^Vw{l) for large I, we can make 
use of a renormalization-group arguments to determine the allowed behaviour. As h 
and tg (X 9 — a are always relevant operators, we will restrict ourselves to the filling 
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transition boundary 9 — a, h — 0. The effective 2D wetting binding potential decays 
faster than 1 /rj'^ , so it becomes irrelevant in the renormalization-group sense [361 137) . 
The renormalization-group flows are dominated by the unstable fixed point and the 
stable fixed point for the 2D wetting binding potential —5/ 7277^, which correspond 
to the tricritical and critical filling transition, respectively. Thus, for large scales the 
only effect of such binding potentials that decay faster than l/l^ is to renormalize the 
line tension associated with the wedge bottom. 

For long-ranged binding potentials i.e. those for which, at = a, l^Vw{l) — > oo 
for large we may resort to a mean- field analysis [SJIS]. We expect that close to the 
filling transition the PDF is asymptotically a Gaussian characterized by the mean- 
field interfacial height lyy^ and roughness ■ However, if the next-to- leading order 
to the wedge binding potential is determined by the short-distance repulsive part of 
WTr{z), we find a similar scenario to the one depicted in figure [5l In particular, it is 
possible to bind the interface to the wedge bottom. The interfacial roughness will be 
controlled by the (microscopic) lengthscale I* ~ /tt corresponding to the maximum of 
the total effective binding potential. The existence of such a bound state, as well as 
the threshold to the critical filling transition will depend on the specific details of the 
interfacial binding potential. 

The predicted filling phenomenology presented in this paper can hopefully be 
checked experimentally or by computer simulations of more microscopic models. Of 
course our predictions for contact (strictly short-ranged) forces requires the elimination 
of van der Waals forces which are ubiquitous for simple fluids. Here our results are 
most easily tested using large scale Ising model simulations. In this case it should 
be straightforward to induce first-order filling by weakening the local spin-substrate 
interaction near the wedge bottom. In contrast our predictions for first-order filling 
with van der Waals forces may well be amenable to experimental verification sometime 
in the near future. Taking an even broader perspective it may be that the chemical 
decoration of a wedge bottom will provide a practical means of eliminating large scale 
interfacial fiuctuations. This may be of relevance to the construction of microfiuidic 
devices whose efficiency will depend crucially on the control of fiuctuation effects. 
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Appendix A. Evaluation of Z{lb,la\Y) at 6 = a for short-ranged forces. 

The evaluation of the interfacial properties at the filling boundary 9 — a in based 
on the knowledge of the partition function ([6]). In this Appendix we will evaluate 
explicitly the partition function Z{lb,la,Y) at the tricritical and critical points for 
short-ranged forces, and we will discuss some properties of the partition function for 
arbitrary u. 

Our starting point is the solution of (fT5|) . for A{1) = 2T,l/a and Vw{l) — 0, so the 
effective potential is V^{ri) ~ —h/12rf'. This Hamiltonian has (1,1) deficiency indexes 
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on the (0, oo) interval, so the solution is not unambiguosly defined, but instead depends 
on a parameter c which defines the short-distance behaviour of the eigenfunctions, i.e. 



cv + r(-l/3) ^ 



7] ^ 



We define the partition function Z{rji,^ rja, Y) via the spectral expansion: 



(A.l) 



(A.2) 



where the summation must be understood as an integral for the continuous part of 
the spectrum. From (I17p . this partition function is related to Z{lb,la,Y) via: 



/ 25] 

ZikJa,Y) = J — {lbl,y/*Zirjik),v{la),Y) 
V a 



(A.3) 



The continuous part of the spectrum of any self-adjoint extension of our interfacial 
Hamiltonian corresponds to £^ > 0. The corresponding scattering states can be 
expressed as [55] : 



Vc2-c(2£;)i/3 + (2£;)2/3 

where Js{x) is the s-order Bessel function of first kind. 
In addition, for c > there is a bounded eigenstate: 



33/2^^3^ 



(A.4) 



(A.5) 



with associated eigenvalue Eq = where Ks{x) is the s-order modified 

Bessel function of second kind. Transforming back to the original variable / = 
(9a/8E)^/^77^/^, the associated eigenfunction ipoil) is: 




ipo{l) cx 



where Ai(a;) is the Airy function, and f„ is defined: 



/ a \ i/3 1 
V2Ey c 



AY.E0 



-1/3 



(A.6) 



(A.7) 



On the other hand, the longitudinal correlation length along the wedge axis reads 
= 2/c^. Comparison with the results obtained by the regularization procedures in 
the text allows us to identify c cx (u — Wc). 

For general c, we cannot perform the spectral integral. However, for c = 
(u = Uc) and c = ~oo (u <C Uc) we have closed form expressions for Z[rib,ria,Y) 

Eaiioi: 



dEe 



Y 



exp 



J_i{V2Ef]a)J^i{V2Er]b) 
1 



I 



(^) 



Zc=^oc — y/VaVb 



2Y 

dEe''^^ Ji{V2Er]a)JiiV2ETjb) 



y/VaVb 

Y 



exp 



2Y 



(^) 



(A.; 



(A.9) 
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with the scaling property: 

Zc=o,-o.^^Uc(^,^] (A.IO) 



/Y \VY VY 
Substituting into (jA.3|) we obtain 

~ 36^^ I — W-* — ) ^-i [-^Y^ ) ^^-^^^ 

where I* = l/£_u and Y* = Y/^y. Now the lengthscale ^„ is arbitrary but £,y — 4E^2/"^- 
In both cases we have scaling such that 

for c — and —oo. For Y — > oo, Uc has the following asymptotic behaviour: 



Uc^o (A.14) 

7 \ 3/2 / , \ 3/2 



oo « (^^j (^^j (A.15) 

To obtain results for arbitrary c, we will make use of the Krein formula ^41^ , which 
relates the Green functions of different self-adjoint extensions of a closed symmetric 
operator. In our case, the Green function Z{lb, la] E) is basically the Laplace transform 
of 2'(/f,, la \ Y) with respect to Y: 

/■oo 

Z{lb,la;E)^ / dYexp{EY)Z{lbJa;Y) (A.16) 
Jo 

For c — and c = — oo we have the closed expressions for the Green function [39j : 

2c=o = ^-^K. ^y-E*{l*^r^ I_. ^y-E*il*^r^ (A.17) 

Zc=^^ = ^-^K. gyr^^) /. g^-i?*G*)3^ (A.18) 
where It, and are the largest and smallest between I* and l^, respectively, and 

E*=E^y. 

To continue, we define the lengthscales and for each c as: 

& = i|j (A-20) 

Note that and £,y reduce to the relevant correlation lengthscales for c > 0. 
Application of the Krein formula leads to the following expression for the Green 
function corresponding to an arbitrary c: 

Zc{lt,la,E)^Z-ooilb,la,E) 



+grjlK. -^-E*[llY\ -J~E*{llf (A.21) 
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where g = g{E, S, a) is obtained by imposing that the short-distance behaviour of Z 
is consistent with the boundary condition (jA.ip . After some algebra, we obtain the 
following expression: 

where the negative sign corresponds to c > 0, and the positive sign to c < 0. It is 
worthwhile to note that the dependence on the boundary condition, i.e. c, has been 
absorbed into the lengthscales ^„ and ^j,. We can formally invert (|A.22p : 

ZcilbJa,Y) = — dEeM~EY)ZSb,la,E) (A.23) 

with 7 < min(0, —c^/2). Note that there is a branch point of at £' = 0. In addition, 
there is a single pole for E — 
lead to the following scaling: 



V2 for c > 0. The scahng properties of Zc (|A.22p 



(A.24) 



(jA.23[) provides a means of obtaining the asymptotic behaviour oi Zc a,s Y —> oo. 
For c > 0, the dominant contribution comes from the pole of Zc a,t E = —c^/2, so: 



Zcilb,la,Y) 



/7I7*Ai(/:)Ai(/:)exp(r*) 



(A.25) 



in complete agreement with our previous analysis. Thus and are true correlation 
lengths. 

For c < 0, the branch point of Zc at i? = controls the large-l" behaviour of Zc- 
We expand Zc around E = 0: 




(A.26) 



Thus Zc has for large Y the following asymptotic behaviour: 



7* /* 



e„r(-i) (r*)4/3 
3^^^r (1) 




32/3r(i 



+ .. 



(A.27) 



Consequently, we find that for lUl > 1 but \/l^/{Y*)^/^ small, the asymptotic 
behaviour of Zc is the same as for the c — — oo case, i.e. critical filling. For the c < 
case £,u and S^y are not correlation lengths, but the (microscopic) lengthscales for / 
and Y, respectively, so the critical filling behaviour is observed over larger scales than 
these. 

In order to obtain the asymptotics of the correlation function h{lb, la, Y) for c > 0, 
we also need the next-to- leading order dependence of F that emerges from Zc{lb, la,Y). 
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Our approach provides a systematic means of obtaining such correction which can be 
obtained in a similar manner as before. If we denote by the leading contribution 
to Zc given by (jA.25[) . the first correction is given by: 
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